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Abstract 



We investigate the Big-Bang Nucleosynthesis constraints on neutrino mixing 
in the framework of the two four-neutrino schemes that are favored by the re- 
sults of neutrino oscillation experiments. We discuss the implications of these 
constraints for terrestrial short and long-baseline neutrino oscillation experi- 
ments and we present some possibilities of testing them in these experiments. 
In particular, we show that from the Big-Bang Nucleosynthesis constraints 
it follows that the k^— > i> T transition is severely suppressed in short-baseline 
experiments, whereas its oscillation amplitude in long-baseline experiments 
is of order 1. We also propose a new parameterization of the four-neutrino 
mixing matrix U which is appropriate for the schemes under consideration. 

14.60.Pq, 14.60.St, 26.35.+C 
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I. INTRODUCTION 



The problem of the masses and mixing of neutrinos (see Refs. ]I]-0) is the first pri- 
ority problem of neutrino physics. Many experiments searching for neutrino oscillations 
and neutrinoless double-beta decay and investigating the high-energy part of the tritium 
beta-spectrum are going on or are under preparation. At present, indications that neutri- 
nos are massive and mixed have been found in solar neutrino experiments (Homestake ||, 
Kamiokande |J, GALLEX [[7], SAGE [|] and Super-Kamiokande |],[nj), in atmospheric neu- 



trino experiments (Kamiokande |TT|, 1MB ||12|| , Soudan |13|| and Super-Kamiokande flUyiOl ) 



and in the LSND experiment |]T5|. From the analyses of the data of these experiments 
in terms of neutrino oscillations it follows that there are three different scales of neutrino 
mass-squared differences: 

Am s 2 un ~ KT 5 eV 2 (MSW) or Am s 2 un ~ KT 10 eV 2 (vac. osc.) [16,0 , (1.1) 

Am 2 tm ~5xl(T 3 eV 2 @ , (1.2) 



Am 2 BL ~leV 2 0, (1.3) 



where Amg BL is the neutrino mass-squared difference relevant for short-baseline (SBL) ex- 
periments, whose allowed range is determined by the positive results of the LSND experi- 
ment. The two possibilities for Am 2 un correspond, respectively, to the MSW |L9| and to the 
vacuum oscillation solutions of the solar neutrino problem. 

Thus, at least four light neutrinos with definite masses (with only three active flavor 
neutrinos, u e , z/ M , v T ) must exist in nature in order to accommodate the results of all neutrino 
oscillation experiments. This means that there exists at least one non-interacting sterile 
neutrino [^Oj -pBj- m this case, a left-handed neutrino field u aL (a = e,/i, r, s) is a mixture 



of the left-handed components v\a, of the four fields of neutrinos with definite masses rrik 
(* = 1,...,4): 

4 

VgL = Uak VkL (tt = 6, /i, T, s) , (1.4) 

fe=l 

where U is the 4x4 unitary mixing matrix. 

In Ref. |22[] we considered the schemes with four massive neutrinos and we have shown 
that from six possible mass spectra of four massive neutrinos only the following two mass 
spectra with two pairs of close masses separated by a gap of about 1 eV (the "LSND gap") 
are compatible with all existing data (see also Ref. p3|): 



atm sun sun atm 

(A) mi < mi TOs < ml and (B) mi < ml < ml . (1.5) 
v ^ , v ^ / 

SBL SBL 

In scheme A, Am 2 tm = m\ — m\ and Am 2 un = m\ — m\ are the mass-squared differences 
relevant for atmospheric and solar neutrino oscillations, respectively, whereas in scheme B 
we have Am 2 un = m\ — m\ and Am 2 tm = m\ — m\. In both schemes Amg BL = m\ — m\ is 
relevant for neutrino oscillations in short-baseline experiments. 
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In this paper we study the compatibility of these two schemes with an upper boundQ 
N v < 3.9 (see Ref. f27j) on the effective number N v of light neutrinos relevant in Big-Bang 
Nucleosynthesis (BBN) (see, for example, Ref. |28[). We derive the constraints on some 
elements of the neutrino mixing matrix U that follow from N u < 3.9 (see also Ref. f23|) and 
we discuss some possibilities to check these constraints in future short-baseline (SBL) and 
long-baseline (LBL) neutrino oscillation experiments. 

Let us notice that there was recently some discussion in the literature about the validity 
of the upper bound N v < 3.9 (see Ref. f27|). This discussion was initiated by conflicting 
results of different measurements of the deuterium abundance in high-redshift hydrogen 



clouds p9|j30|| . The situation is not completely clarified, but there seems to exist now a 
strong case in favor of N v < 3.9 |27[ . 
Let us define the quantities 



(a = e, /x, r, s) . (1.6) 



fc=l,2 



Okada and Yasuda [2]J have shown that the bound N v < 3.9 implies stringent limits on c s in 
scheme A and 1 — c s in scheme B. This means that the solar neutrino deficit is explained by 
small angle MSW v e — > v s transitions^ and the atmospheric neutrino anomaly by — > v T 
transitionsf] (see Ref. ||18||). In this paper we have reanalyzed the constraints on the parameter 
c s that follow from BBN. We have obtained more stringent upper bounds on c s in scheme A 



and 1 — c s in scheme B than those presented in Ref. PBI and we have investigated the effects 



of these constraints for SBL and LBL neutrino oscillations, some of which can be tested in 
the near future. 

The plan of the paper is as follows. In Section |J we summarize the formalism of SBL 
neutrino oscillations in the schemes A and B (for more details see Refs. [T2T],£?2|). In Section 
III] we investigate in detail the dependence of N v on c s . From the BBN upper bound on N„ 
we obtain upper bounds on c s in scheme A and 1 — c s in scheme B. In Sections fV| and |V| we 
obtain various general relations between SBL and LBL oscillation amplitudes and c s . These 
relations can be useful to check the BBN constraints and to test the hypothesis of existence 
of a sterile neutrino. In Section [V] we also introduce a new parameterization of the mixing 
matrix that is suitable for the schemes under consideration. 



1 For N u > 4 the schemes A and B are trivially compatible with BBN. 

2 The vacuum oscillation and the large angle MSW solutions of the solar neutrino problem are 
disfavored by the data in the case of v e — > v s transitions (see Refs. [3 1,3^ Jill). 

3 This explanation of the atmospheric neutrino anomaly is consistent with the recent results of 
the first long-baseline reactor experiment CHOOZ |33| in which no indications in favor of v e dis- 
appearance were found. 
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II. NEUTRINO OSCILLATIONS IN SBL EXPERIMENTS 



In the schemes under consideration the probability of v a — > v@ transitions {(3 7^ a) and 
the survival probability of v a in SBL experiments are given by f21 



1 



Am 2 SBL L 



P» a ^u = 2^«;/3 I 1-cos 

1 / AmH RT L 
Pu -+ v = 1 - - B a . a 1 - cos 

2 ' V 2 p 



with the oscillation amplitudes 



(2.1) 
(2.2) 



a;/3 



B r 



E u * u « 

i 




5>. 



where the index i runs over the values 1, 2 or 3, 4. 

In the following, we will consider values of Awig BL in the wide range 

lO^eV 2 < Awig BL < 10 3 eV 2 . 



(2-3) 



(2.4) 



(2.5) 



The direct comparison of the regions in the A^ e -/S>m\ Blj plane allowed at 90% CL by the 
results of the LSND experiments with the regions excluded at 90% CL by the results of 
other SBL experiments indicate a considerably smaller allowed range for Amg BL : 



0.27 eV 2 < Am 2 T <2.2eV 2 . 



(2.6) 



The results of the combined analysis of the data of all SBL experiments presented in Ref. |3 
show that the lower bound Amg BL > 0.27 eV 2 is robust (if the LSND indication is correct), 
whereas values of Amg BL larger than 2.2 eV 2 maybe are not excluded. 

Taking into account the results of solar and atmospheric neutrino experiments, from the 
exclusion curves of SBL disappearance experiments at each possible value of Amg BL we have 



22 



and 



c e < di 



and 



in scheme A 
in scheme B . 



where 



5" 



[a 



e,/i) 



(2.7) 
(2.8) 



(2.9) 



and B®. a is the upper bound for the amplitude B a;a that can be obtained from the exclusion 
plots of SBL reactor and accelerator disappearance experiments. From the 90% CL exclusion 

v e — > v e reactor experiment and of the CDHS 



curves of the Bugey |35 



61 and CCFR 37 



v ii ~^ v \i accelerator experiments it follows that a° e is small {pP e < 4 x 10 2 ) 
the range ( p~5|) and a? < 10 _1 for Amg BL > 0.5 eV 2 (see Fig. 1 of Ref. 



for Amg BL in 
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III. BBN CONSTRAINTS ON 4-NEUTRINO MIXING 



In this section we derive the BBN constraints on the elements of the mixing matrix in 



the two schemes A and B. We follow the standard arguments p9|- ^2| , p3| leading to the BBN 
constraints on the elements of the mixing matrix U if a sterile neutrino exists. According to 
the standard BBN scenario (see, for example, Ref. |p8| ) a lepton asymmetry can be neglected 
and the active neutrinos are in chemical equilibrium. Thus, the following considerations 
are valid for a temperature range T dec < T < 100 MeV with the decoupling temperature 
^dec ~ 2 MeV for electron neutrinos and T& ec ~ 4 MeV for muon and tau neutrinos. In this 
temperature region, the effective potentials of neutrinos due to coherent forward scattering 
in the primordial plasma are given by |43|j4l|j4! 



V e = V = -Q.02G F p 



rp4 



V^r = £V and V 8 = 



(3.1) 



where p is the neutrino momentum, which henceforth will be replaced by its temperature 
average (p) ~ 3.15 T, and £ = cos 2 6^/(2 + cos 2 6^) ~ 0.28. The propagation of neutrinosfj 



is governed by the effective hamiltonian [IE 

1 



2p 



U diag(m 2 , ml m 2 , m 2 ) U~ l + diag(^, ^V, £V, 0) 



U'di ag (E 1 ,E 2 ,E 3 ,E 4 )U 



i-i 



(3.2) 



where Ex, E 2 , E 3 , E^ are the effective energy eigenvalues of neutrinos in matter and U' is 
the corresponding effective mixing matrix. The amount of sterile neutrinos present at nu- 



cleosynthesis can be calculated using the differential equation p9 | 

q ^ ] {Pva^Vs) C0\{^ V a (1 _ W V S ) 



dt 



(3.3) 



where n Va is the number density of the sterile neutrino relative to the number density of 
an active neutrino in equilibrium and T Va are the collision rates of the active neutrinos, 
including elastic and inelastic scattering | 4"2"|| , given by 



4.0 G 2 F T 5 and T, 



BT U with 5 = 0.722. 



(3.4) 



Let us notice some important characteristics of Eq.(p.3|): 



i. The probabilities for v a — > v s transitions are averaged over the collision time t co \\ 
1 /r„ e and also n Ug has to be considered as such an averaged quantity. 



4 In the absence of a lepton asymmetry the neutrinos and antineutrinos evolve identically and 
hence only neutrinos are considered here. 
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ii. This equation should approximately hold for non-resonant and adiabatic resonant 
neutrino transitions, as long as the active neutrinos are in chemical equilibrium and 
W "C tcoii <S t exp is fulfilled. The characteristic expansion time of the universe t exp 
is given by t exp = 1/H where H is the Hubble parameter, which is related to the 
temperature T by H = -f/T ~ 0.69 (T/l MeV) 2 s -1 . 

iii. The condition t osc <C t co \\ means that neutrino oscillations have to be fast relative to 
the collision time. The relation T Ue /H ~ 1.2 (T/l MeV) 3 shows that for temperatures 
larger than 2 MeV the collision time is always much smaller than the expansion time 



iv. Chemical equilibrium for the active neutrinos is maintained at temperatures larger 
than the neutrino decoupling temperature even in the presence of sterile neutrinos, 
provided the relevant mixing angles in the primordial plasma are not too large W2 . 



Using the relation H = —T/T, the evolution equation ( |3.3j ) for n Us can be rewritten as 
dn,. 



dT 2HT 



\?p Yl ( P ^^)coiir^(l -n Us ) . (3.5) 



Since by definition N u is the effective number of massless neutrino species at Td ec , in order 
to get a constraint on the mixing of sterile neutrinos we need to calculate the value of n Ug 
at T doc produced by neutrino oscillations. With the initial condition n Va (Ti) = (T ~ 100 
MeV), the integration of Eq.(|3~5|) gives [HI 



tU (Pv^» s ) con? v Jt) =l- eX p(-F). (3.6) 



T, 




Imposing the upper bound 

n„ a (T dec ) <SN = N u -3, (3.7) 

one obtains the condition 

F < | m(l — 5N) \ . (3.8) 



Since c e is small in scheme A and 1 — c e is small in scheme B (see Eqs. fl2.7D and 
they do not have any effect on neutrino oscillations during BBN. Hence in this Section we 
use the approximation c e = in scheme A and the approximation 1 — c e = in scheme B 
23| . In this case, the elements U e i and U S i (i = 1, . . . , 4) of the mixing matrix can be chosen 



real and can be parameterized by 

(17^) = (0, 0, c e , sg), (U si ) = (s^Sx, -s^Cx, -c^se, c^cq) in scheme A, (3.9) 
(U^) = (ce, se, 0, 0), (U si ) = (s^sg, -s v c e , -c v s x , c^c x ) in scheme B, (3.10) 

with the abbreviations c p = cos p and s p = sin p for p = 9,x,f and < ip < it/2. In both 
schemes A and B we have c s = sin 2 <p. 
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Then, it can be shown that 

/ j V-fl/a-^I/s/COll-L ^ a — 2 



1 - 2D£x atm cos2x + ^fetm) 2 



4c s (l — c s ) . 
+ 1-2^ S bl(1-2c s ) + (^sbl) 2 ' ' 1 ' ! 



where 

^Sr- 218 * 10 "^ (l4vj (-atm.SBL) (3.12) 
and p is replaced by its temperature average. With the definitions 

D = c s in scheme A and D = 1 — c s in scheme B , (3.13) 



Eq. ( |3.11| ) holds in both schemes. The two terms in Eq. ( |3.11|) correspond to oscillations due 



to Am^ tm and Amg BL , respectively. Oscillations due to Am^ are neg lected in Eq. (|OlD 



because their contribution to F (see Eq.(|3.6|)) is at least two orders of magnitude smaller 
for T > T dec . 

Using the expression ( |3.11| ) one can find an analytic form of the integral ( |3.6| ) and the 



condition (|3.8|) gives the bound 

920 ( Cs V /T^ + 33 ( An foO V2 - Sin2 2x D" 2 < | ln(l - 8N) \ . (3.14) 
V leV 2 J \10- 2 eV 2 J Vl + cos2x ~' V " V ; 

With the definitions ( |3.13j ) this bound holds for both schemes, provided the conditions laid 
down in i.-iv. are fulfilled. 

Let us now discuss the implications of the bound ( |3 . 1 4| ) . In order to understand if this 
bound is compatible with a large D, we first consider the possibility that D is close to one. 
In this case we have {U^i.U^) ~ (cos x, sin x) in scheme A and {U^,U^) ~ (cos x, sin x) 
in scheme B. This means that, in order to accommodate the atmospheric neutrino anomaly, 
sin 2 2% cannot be small. However, this is in contradiction with the inequality ( 3.14 ) because 
of the second term in the left-hand side. Therefore, we conclude that the bound ( |3.14| ) 
implies that D is small, i.e. c s is small in scheme A and 1 — c s is small in scheme B. 

Let us now check the validity of this reasoning considering the oscillations due to Am 2 tm 
which generate the second term in Eq. (|3.14j ) and taking into account that the method 
considered here only holds for non-resonant neutrino oscillations or adiabatic transitions 



through a resonance. Looking at Eq. ( |3.11|) we see that there is a resonant behaviour of 
the oscillations due to Am 2 tm at the temperature where the relation Dt;x a tm = cos 2% is 
satisfied. Since x atm is negative, this is only possible if cos 2% is negative. However, since we 
do not know the sign of cos 2%, we have to check if the resonance is crossed adiabatically in 
the case of a negative sign. A measure of the adiabaticity of the resonance is given by the 
parameter 



<?atm 



A-E atm dt 



2 \ -1/2 
atm | 

10- VV 2 ; v^D(l + cos2x) 3 / 2 



1.23X10-M Am ^X = Sin2 * (3.15) 
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where AE atm and x' are the oscillation frequency and the effective mixing angle in matter. 
If Qatm ^ 1 the adiabatic crossing of the resonance is guaranteed. Clearly, for D ~ 1 and 
sin 2 2x ~ 1 this condition is fulfilled. Furthermore, oscillations due to Am 2 tm occur around 
the temperature T atm ~ ID^ 1 ^ MeV, where D|£x atm | ~ 1. One can show that for D ~ 1 
the relation t osc <C t co n holds for resonant and non-resonant oscillations due to Am 2 tm and 
thus the conditions for the validity of Eq. ( |3.3|) and all relations derived from it are fulfilled 
as well. Consequently, the exclusion of the case D ~ 1 with the help of the second term 
in Eq.( |3.14D is correct. It is interesting to note that no conclusions can be drawn from the 
second term in Eq. ( |3.14j ) if D < 0.1 because then T atm > 10 MeV and in this temperature 
region one has t osc 3> t co n and therefore no oscillations due to Am 2 tm are possible because 
of strong quantum damping. 

Let us now examine the implications of the first term in Eq. ( |3.14j ). The presence of this 
term implies that either c s or y/1 — c s has to be very small. As we have just shown considering 
the second term in Eq. ( |3.14 ), in scheme A c s is small. This implies that 1 — 2c s > and 
only non- resonant oscillations are possible (see the last term in Eq. (|3.11|) ). Hence, Eq.( p.l4|) 
is valid in scheme A and the first term gives the bound 



c, < 1.1 x 10" 



/ Am 2 



SBL 



V leV 2 



-1/2 



ln(l-5iV)| = B W . 



(3.16) 



On the other hand, considering the second term in Eq. (|3.14|) , we have shown that in 
scheme B the quantity 1 — c s is small. Therefore, in scheme B we have 1 — 2c s < and 
a resonance occurs at the temperature T res = 16(Am| BL /l eV 2 ) 1/6 |l - 2c s | 1/6 MeV, where 
^sbl = 1 — 2c s . One can show that this resonance is not passed adiabatically, hence the 
conditions for the validity of Eq. (p.3|) are not fulfilled and the first term of Eq.( |3.14| ) does 
not apply. In this case the amount of sterile neutrinos produced at the resonance through 
non- adiabatic transitions can be calculated with the formula [55,56] 



1 
2 



^ - Plz ) <*os2 r Yros2,;\, ~ .1 



LZ 



(3.17) 



where iff, (ip a ) is the effective mixing angle before (after) the resonance. The last approxima- 
tion in Eq. (|3.17|) follows from the fact that cos 2(pb ~ 1 because of the high effective neutrino 
potential before the resonance and cos2<^ a ~ cos2<£> = 1 — 2c s ~ — 1 for small 1 — c s . The 
quantity Plz = exp(— Q) is the Landau-Zener probability [^Tj and Q is given by [[ij,[4j,[4j 



2tt c s (1 - c s 
h II -2c, 



y ( rfSSBlA 1 

V dt ) 



9.2 x 10 4 



cJl 



II - 2c s | 3 / 2 V leV 



Am SBL 



1/2 



(3.18) 



where St is the half- width of the resonance defined in Ref. [^JJ. Demanding again that 
n Vs < 5N one obtains the condition Q < | ln(l — SN)\ which gives the bound 

-1/2 

1-c, < 1.1 x 10~ 5 ( ^'". SB o L ) |ln(l -SN)\ = B {B) (3.19) 



leV z 
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in scheme B. The linear approximation of the potential used to derive the Landau-Zener 
formula is valid if the potential changes slowly in the resonance region. At the resonance 
one obtains \V/(dV/dt)\/5t ~ |1 — 2c s |/2a/c s (1 — c s ) which is large for small 1 — c s . Hence 
the linear approximation should apply. 

In the following table we show the values of the bounds B^ and B^ obtained from 



Eqs.([3l6|) and fU9p and the LSND lower bound A 



SBL 



> 



0.27 eV 2 (see Eq.Q) for 



different values of 5N: 



SN 


0.9 


0.8 


0.7 


0.6 


0.5 


0.4 


0.3 


0.2 


scheme A: 


B( A )/10- 3 


4.8 


3.4 


2.5 


1.9 


1.4 


1.1 


0.7 


0.5 


scheme B: 


B^/10- 5 


4.9 


3.4 


2.6 


1.9 


1.5 


1.1 


0.8 


0.5 



(3.20) 



The bounds ( |3.16 ) for the non-resonant oscillations in scheme A and ( |3.19 ) for the resonant 
case in scheme B are in rough agreement with the corresponding numerical integrations 
of the evolution equations of an ensemble of oscillating neutrinos which take into account 



2-neutrino oscillations [|4*2||[6" |. 

Let us now summarize our findings. The constraints on the parameter c s in scheme A 
(Eq. ( |3.16| )) and 1 — c s in scheme B (Eq. fl3.19"D ) has been derived in both schemes from the 
oscillations due to Amg BL at temperatures of the order T ~ 16 MeV where £|xsbl| ~ 1- 



The second term in Eq. fl3.14|) , which is generated by the oscillations due to Am^ tm , only 
serves to exclude large D. The transitions of active into sterile neutrinos due to Awig BL are 
non-resonant in scheme A, whereas they are resonant in scheme B, leading to a stronger 
constraint on 1 — c s than the constraint on c s in scheme A (see table ( 3.20| )). We have 
obtained stronger bounds than Ref. [^3| because we have taken into account the complete 
collision rates 



presented in Ref. ||42|| , which are nearly an order of magnitude larger 
than those used in Ref. |23j|. Furthermore, our method shows that the dependence of the 



bound on c s in scheme A and 1 — c s in scheme B on the allowed number of neutrino species 



N u is logarithmic |44[ . 

We want to emphasize that our calculation of the amount of sterile neutrinos brought 
into equilibrium by oscillations before the onset of BBN is based upon equation ( |3.3|) in 
scheme A and upon the Landau-Zener approximation in scheme B. Therefore, the numbers 
given in table ( |3.20| ) should be viewed as order of magnitude estimates of the upper bounds 
on the parameter c s in scheme A and 1 — c s in scheme B. 

Concluding this section, we would like to make some comments on the non-standard 
BBN scenario with a non-zero lepton asymmetry presented in Ref. ||48| . There it has been 



shown that for a certain range of the mixing parameters and for resonant transitions due to 
Amg BL a lepton asymmetry as small as 10~ 10 can be amplified up to 10 -2 and in this way 
transitions at a lower temperature (like the transitions due to Am 2 tm in the schemes under 
consideration) can be prevented. Such a scenario would not change our bound for 1 — c s in 
scheme B, but could have some effect on our considerations for scheme A. It c s is large in 
scheme A, a resonance due to Amg BL occurs and, if a large lepton asymmetry is generated, 
the oscillations due to Am 2 tm are suppressed and the second term on the left-hand side of 



Eq. fl3.14|) is not valid. Hence a large parameter c s cannot be excluded. Therefore, contrary 
to our discussion in this section, in this case we would have two possibilities for scheme A: 
either c s is small and has to obey the bound ( 3. 16] ) or c s is close to one and 1 — c s has to 
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obey the same bound as given in Eq. (|3.19|) for scheme B. However, judging from the results 
of Ref. it has yet to be shown that with c s close to one a realistic scenario meeting all 
experimental constraints can actually be achieved in scheme A and we leave this problem 
to future investigation. 

IV. TERRESTRIAL NEUTRINO OSCILLATION EXPERIMENTS AND STERILE 

NEUTRINO MIXING 

As we have shown in the previous Section, BBN suggests that the parameter c s in scheme 
A (1 — c s in scheme B) is small. In this Section we present some possibilities to obtain 
information on the parameter c s from future results of SBL and LBL neutrino oscillation 
experiments. We consider explicitly scheme A, but the same results are valid in scheme B 
after the replacement of c s with 1 — c s . 

Let us start with the consideration of — > v T oscillations, that are presently searched 
for in the CHORUS [3^] and NOMAD |50| experiments and will be searched for in the 
COSMOS |H] experiment. From Eq.( |2.3| ), for the amplitude of — ■> v T oscillations we have 
the upper bound 

A KT <4(l-c M )(l-c r ). (4.1) 
The unitarity of the mixing matrix implies that 

J2 c « = 2 ( 4 - 2 ) 



and Eq. ( (4.1[ ) can be rewritten as 

Afi-T < 4(1 - c M ) [ce + c,-(l- c„)] . (4.3) 
Since the quantity 1 — satisfies the inequality ( p.7|) , if 

c e + c s > 2al (4.4) 

we have 

A„., T < Aa%c e + c s - a°) < 4aJ(a° + c s ) . (4.5) 
On the other hand, if the inequality 

c e + c s < 2aJ (4.6) 

is satisfied, for the amplitude A^ T we have 

A^ T <{c e + c s ) 2 <A{alf . (4.7) 

As it is seen from Eqs.( |4.5| ) and ( (4.7| ) , if it will be found that the amplitude A^ T satisfies 
the inequality 
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it will mean that the parameter c s must satisfy the inequality 



c s > 2al 



(4.8) 



(4.9) 



and B° (dashed curve 
and CCFR W 



and is in general not small. In Fig.|I] we have plotted the values of 4(a°) 2 (solid curve) 

obtained from the 90% CL exclusion curves of the CDHS |36| 
experiments in the range (|2.5| ) of Amg BL , that covers the LSND- allowed 
region. The dashed curve representing B ^ constitutes an upper bound for A^. T due to 
the conservation of probability. Therefore, if the inequality ( |4.8|) is satisfied, A^ T must lie 
between the solid and dashed curves. In Fig.^we have also plotted the most recent exclusion 
curve presented by the CHORUS collaboration |52| (dash-dotted curve), the expected final 
sensitivity of the CHORUS j|9[ and NOMAD j50| experiments (dash-dot-dotted curve) and 
the expected sensitivity of the COSMOS experiment (dotted curve). One can see that 
a substantial part of the region in which the inequality (|4.8|) is satisfied is already ruled 



out by the present CHORUS exclusion curve, a large part will be excluded when the final 
sensitivities of the CHORUS and NOMAD experiments will be reached and almost all the 
region will be excluded if the COSMOS experiment will not find — > v T oscillations. 

t2, obtained from the 90% CL exclusion 



In Fig.|2| we have plotted the values of 2a° 
curves of the Bugey [0, CDHS M and CCFR 



values of Amg BL in which a° e < 2a°, if the inequality 



experiments. One can see that for the 
( ^4.8| ) is satisfied c s must be large and 
incompatible with the BBN upper bound. On the other hand, for the values of Awig BL in 
which a° e > 2a°, c s does not need to be large (the inequality ( f4.4| ) is satisfied because c e can 
be larger than 2a°) and there is no contradiction with the BBN upper bound on c s even if 
the inequality ( |4.8| ) is satisfied. If A^. T will be measured and found to satisfy the inequality 
( |1.8| ), also Amg BL will be determined with some accuracy. In this case, depending on the 
allowed range of Amg BL , it will be possible to decide if a small c s compatible with the BBN 
upper bound is excluded or not. Note, however, that for the LSND range (|2.6|) of Am| BL 
the parameter c s must be large if the inequality ( f4.8p is satisfied. 

Furthermore, if some future experiment will measure a value of A^. T such that Eq.( 



is satisfied, then the inequality ( |4.5|) will imply in addition the lower bound 



c s > 



A 



ff.T 



- a: 



(4.10) 



if A KT > 



,0\2 



which is more stringent than the lower bound (|4 

Let us now consider the possibility that a future experiment will measure a small value 
for the amplitude A^. T , i.e. such that the inequality ( |4.7|) is satisfied. This measurement will 
not give any information on the value of c s . In this case, the BBN bound derived in Section 
|H implies that c s ~ and for the amplitude A^ T we have the upper bound 



A,..- < (a 



0\2 



(4.11) 



The numerical values of this upper bound in the range ( |2.5| ) of Amg BL obtained from the 
90% CL exclusion curve of the Bugey [^] experiment is shown in Fig.[| (solid curve), to- 
gether with the most recent exclusion curve presented by the CHORUS collaboration 
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(dash-dotted curve), the expected final sensitivity of the CHORUS |49[ and NOMAD |50 
experiments (dash-dot-dotted curve) and the expected sensitivity of the COSMOS |[I ex- 
periment (dotted curve). 

Thus, we have seen that the investigation of — > v T oscillations in the region of Amg BL 
that includes the LSND-allowed region could allow to obtain information on the value of the 
parameter c s that is important for BBN. 

Let us consider now the possibility to obtain information on the parameter c s from the 
data of inclusive — > and v e — > v e experiments. From Eqs. (|2.4j ) and ( |2.7| ) it follows that 

1 



Ce = -[l-y/l-B e]e ) , (4.12) 
1 



1-S= 2 t, 1 - V1-^;mJ • (4-13) 
Furthermore, from Eq. (|4.2j ) we have 

c s >(l-c^)-c e . (4.14) 
From Eqs. ( |4.12|) — (|4.14| ) we obtain the inequality 

c s > \ (a/1 " B e]e - y/1 - B^j . (4.15) 

The amplitude B e . e is small in the whole range Q2.5|) of Am| BL . If the parameter c s is small 
as suggested by BBN, also the amplitude must be small. The existing data do not 
exclude, however, large values of B^ for Am| BL < 0.3 eV 2 . Hence, our analysis shows that 
the investigation of inclusive — * transitions for Amg BL < 0.3 eV 2 would be interesting 
for the check of the constraint that follows from BBN. 

Let us now consider the amplitude A e - S of v e — > v s oscillations. From the upper bound 

A e]S <Ac e c s <4a°c s , (4.16) 

which follows from Eqs.( |1.6|) , (|2.3| ) and ( p.7|) , it is clear that the BBN bound on c s implies 
that A e - S is extremely small. In this case the amplitudes measured in SBL v e disappearance 
experiments and SBL v e —* and v e — > v T appearance experiments are related by 

B e . e ~A me + A e . T . (4.17) 

The experimental check of this relation allows to test the BBN bound on c s . Indeed, if 
SBL v e disappearance experiments will find neutrino oscillations with a lower limit B^ a ' 
for the amplitude B e . e such that -E>i™ m ' ) > A^e + Ai™ , where A^e and Ae™ ax ' ) are the 
experimental upper limit for A^ e and A e . T , respectively, it will mean that A e . s > 0. In this 
case, taking into account that B e - e = A^ e + A e - T + A e . s and that for small B e - e we have 
c e ~ B e . e /4:, the inequality A e - S <4c e c s gives the lower bound 

A A 4- A /|( max ) I 4 (max) 

c > 251 ~ i _ A * e + Ae ' r > i _ +/le;r (a is) 

s ~ 4 c B ~ R ( min ) ' 1 j 
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If, for example, it will be found thatQ B { ™ n) ~ 10" 2 , A^ ax) ~5x 10~ 3 and A^T^A^, 
Eg. ( |4. 18| ) will imply the lower bound c s > 0.5, which is incompatible with the BBN upper 
bound. 

CP violation in long-baseline neutrino oscillation experiments can also be used to 
check the BBN constraint on c s , though in practice this possibility is more remote. 
It has been shown in Ref. |25| that the absolute value of the CP-violating parameter 
V = 4Im[C/ /Lt i U* x U* 2 U t2 ] (in scheme A, whereas V = 4 Im [U^ U* 3 £/"* J7 t4 ] in scheme 
B) could be as large as 2/3 \/3, which is the maximal value allowed by the unitarity of the 
mixing matrix. With the method discussed in Ref. one can show that 

| VI < (Ce + C 8 )y/1-C e -C 8 (4.19) 

for small c e + c s . From the upper bound c e < aP e (see Eq. (|2.7f )) and the BBN upper bound 
c s < B^ (see Eq.( |3.16| )) we have 



I VI < (a° + B (A) )y 1 - a° - B^ . (4.20) 

Consequently, taking into account that a° is small, the BBN upper bound on c s implies that 
CP violation in the — > v T channel is suppressed. On the other hand, if CP violation will 
be observed in LBL — > u T experiments, using Eq. ( |4.19|) it will be possible to set a lower 



bound for c s . Analogous results are valid in scheme B with the replacements c e — > 1 — c e 
c s -> 1-c, and B^ B™. 



V. A PARAMETERIZATION OF THE 4-NEUTRINO MIXING MATRIX 

In this Section we propose a parameterization of the 4x4 mixing matrix U which is 
suitable for the two schemes A and B that allow to accommodate all existing neutrino 
oscillation data. In these schemes the neutrinos v\, v 2 and the neutrinos z/3, 1/4 give separate 
contributions to all observables. Thus, it is natural to consider separately the elements U a i 
with i — 1, 2 and the elements U a k with k = 3,4. In the following we will stick to scheme A 
for the actual presentation. The formulas in scheme B can be obtained with the exchanges 
1 t»3 and 2 4 of the columns of U. 

We take into account that all observable transition probabilities are invariant under the 
phase transformation 

/ — < "" / ' nJ < (5.1) 



5 These values are compatible with the present upper bound for B e - e , which is given by the results 
of the Bugey |35| experiment, and with the present allowed range of A^- e , which is given by the 
results of the LSND |15| and Bugey experiments (B e . e < 4 x 10~ 2 and 2 x 1CT 3 < A^ < 4 x 10~ 2 
for Am| BL in the LSND range fl2.6|) ). Notice that the allowed range of A^- e will be checked in the 
near future by KARMEN |53| and other experiments [|54| . Furthermore, information on A e - T could 
be obtained in the future in experiments using a v e beam from a muon collider (see Refs. [55,56]). 
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where x a and yj are arbitrary parameters. Therefore, U e k with k = 3, 4 and £7" M j with i = 1,2 
can always be considered as real vectors. Let us write them as 



U ek = r e v k _ 2 



(k = 3,4) 

(2 = 1,2), 



with 




U ek \ 2 and 




The unit vectors v and w can be written in the form 

v = (cos#, sin 8) and w = (cos 7, sin 7) . 
Let us introduce also the orthogonal unit vectors 

v — {— sin 8, cos 8) and w 1 = (— sin 7, cos 7) 



(5.2) 
(5.3) 



(5.4) 



(5.5) 



(5.6) 



Then the vectors U a i with a^e and % = 3, 4 can be expanded over the orthonormal basis 
{f , f- 1 } and the vectors U ai with a ^ ji and z = 1, 2 can be expanded over the orthonormal 
basis {w, w^. As a result, we have the parameterization 



U 



( p e w + qeW 1 - r e v \ 
p^w r^v + s^v 1 - 

p T w + q-rW 1 - r T v + SrV 1 - 

\ p s w + qgW 1 - r s v + SgV 1 - J 



(5.7) 



for the mixing matrix. Using the invariance under the transformation (|5.1|) , we can choose 
the parameters q T and s s as real. 

From Eqs. (|1.6|) and (|5.4|) it is obvious that the parameters r e and p M are connected, 
respectively, with the parameters c e and c M by the relations 



Vi - c e , p^ = ye:; . 



(5.8) 



Furthermore, from the unitarity of U it follows that the parameters and p e are related by 

* / c u 

ra = -P e 



1 



(5.9) 



Similarly, we can think of q T and s s being fixed by the unit length of the lines of U and five 
of the remaining complex parameters by the orthogonality of different lines. Counting the 
number of remaining real parameters in U we obtain 10 versus 9 physical real parameters 
in a 4x4 unitary mixing matrix. A careful inspection of the residual phase freedom reveals 
that one of the complex parameters p e , q e , r M , s^, p T , r T , s T , p s , q s or r s can be chosen real 
in addition. 

The significance of the parameters in the matrix ( |5.7| ) shows up by considering the 
amplitudes of neutrino oscillations in disappearance and appearance SBL experiments, which 
are directly connected with the parameters p a and r a . In fact, we have 
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B. 



e;e 



B 



(5.10) 



and 



.4 



We 



|2 2 
I P„> 



4|r, 



|2 2 
\ ' e i 



A 



4|p, 



|2 2 

I ?v 



(5.11) 



The SBL oscillation amplitudes do not contain the angles 9 and 7. This is connected with 
the fact that SBL amplitudes are determined by products of vectors in the 1,2 space or in 
the 3, 4 space (see Eqs.( [2.3|) and ( |2.4|) ). However, information on the angle 7 can be obtained 
from the results of atmospheric and LBL neutrino experiments, whereas the results of solar 
neutrino experiments give information on the angle 9. 

In the four-neutrino scheme A the parameters p e , q e are small because c e = \p e \ 2 + \q e \ 2 



is small (see Eq.i 



Let us now investigate if it is possible to obtain bounds on the 



parameters \p s \ and \q s \ from measurable quantities. We will use the unitarity of the mixing 



matrix and we will work in the approximation c e 
( p.7[ ) of the mixing matrix reduces to 



0. In this case, the parameterization 



U 



/0 v \ 

p^w s^v 1 - 

p T W + qrW 1 - SrV 1 - 



\ p s w + q s w~ 



J 



and the amplitude of — > u T oscillations is given by (see Eq.( |2.3| )) 



(5.12) 



A - 4k l 2 h I 2 



(5.13) 



Furthermore, taking into account that 



I |2 1 I |2 1 I |2 
\Pr\ + \Qt\ + \ S t\ 



I |2 1 I |2 

\Qt\ + \q s \ 



(5.14) 
(5.15) 



we have 



+ \Pr 



2 >\S T \ 2 



(5.16) 



From Eqs. (|5.13|) and ( |5. 16| ) and taking into account that |s M | 2 
we obtain 

1 1 2 \ 

1*1 



This lower bound for \q s 



l-c M <a° (see Eq.Q), 
(5.17) 



= c s — \p s \ 2 , obtained in the approximation aP e = 0, is consistent 
with the lower bound ( 4.10| ) on c s . From the comparison of these two lower bounds it is 
clear that if c s is large, its value approximately coincides with l^l 2 , whereas \p s \ 2 is small 
if a° is small. Indeed, it is possible to obtain an upper bound for the parameter \p s \ 2 using 
the unitarity relation 



PnP*8 + = 0. 



(5.18) 
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Taking into account that 



8^ = 1 -Pt, (5-19) 



s s \ 2 = 1 - \p s \ 2 - \q s \ 2 , (5.20) 



we find from Eq. (|5.18|) 



\Ps 



I \U\- =1 - P »- (5 ' 21) 

Since 1 — p 2 = 1 — c M < a° (see Eq. (|2.7|) ), we obtain the upper bound 

\Ps\ 2 <al- (5.22) 

Therefore, if the experimental upper bound for the amplitude of SBL — > transitions 
is small, |p s | 2 is small as well. 

Note that the observation of SBL inclusive — > transitions together with the obser- 
vation of SBL — > v T and z/ M — > v e (LSND [15|]) transitions would allow to get information 
on the transitions of muon neutrinos into sterile states. In fact, because of the unitarity of 
the mixing matrix, for the amplitude A^ s of — > v s transitions we have 

— B^.p — A^- e — A^- T . (5.23) 

Hence, we think that it is important to continue the search for neutrino oscillations in the 
inclusive — > channel, especially in the region of small Amg BL , below 0.3 eV 2 , which 
has not been explored so far. Let us consider the possibility that the right side of Eq.( |5.23| ) 



will be measured to be different from zero, i.e., that — > v s transitions occur in SBL 
experiments. In the scheme under consideration the amplitude of these transitions is given 
by 

A^ s = Apl\ Ps \ 2 . (5.24) 
Taking into account that p 2 = c M , Eqs.( P~l~3| , ( |5\23D and flOg ) imply that 

i^ffTTfi&r (5 - 25) 

Finally, if it will be found that not only v e — > v e SBL transitions but also — > SBL 
transitions are strongly suppressed, it will mean that the only oscillation channel involving 
sterile neutrinos which may be not suppressed is the channel v T u s . Therefore, in this 
case it will be possible to obtain information on transitions of active neutrinos into sterile 
states only if v T beams (see Ref. [fjTfl ) will be available. In this case |p s | 2 ~ and 

B T , T = A T . s = A\q s \ 2 (l-\q s \ 2 ) . (5.26) 

Let us now consider the possibility that not only p e , q e but also p s , q s are small, as follows 
from the smallness of c s = \p s \ 2 + \q s \ 2 implied by BBN if N v < 3.9, as shown in Section [II. 



In this case it is possible to expand all the quantities p a , q a , r a and s a in powers of the small 
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parameters p e , q e , p s , q s , up to order c e , c s , y/c e c s . This means that in this approximation 
it is possible^ to express all elements of U in terms of the parameters p e , q e , p s , q s , 9 and 
7. Note that one of the phases of the small complex parameters is unphysical and can be 
transformed away reaching thus the number of 9 physically independent real parameters in 
U. Performing this program we obtain 

IpmI 2 - 1 - bel 2 - b s | 2 , r M ~-p*, S„~-p* s , 

Pr ^ -p e q* e - P s q* s , \Qt\ 2 ~ 1 - |g e | 2 - \q s \ 2 i r T ~ -g*, s T ~ -q*, (5.27) 

>"«— ~P*ePs- q* e Qs, |S S | 2 ~1-C S . 

In this approximation, Eqs. (|5.10|) and (|5.11|) become 



B e . e ~ 4(\p e \ 2 + \q e \ 2 ), ~ 4(|p e | 2 + \p s \ 2 ) (5.28) 

and 

A Ke ~4|p e | 2 , A e;T ~4|g e | 2 , ~4|p e g e *+ M *| 2 . (5.29) 
Therefore we have 

B e;e ^ A me + A e . r . (5.30) 
The absolute value of p s is given by 

\Ps\ 2 * \(B^ - A, ;e ) . (5.31) 



In the absence of v T disappearance experiments (see Eq.( |5.26|) ), the determination of the 
parameter |g s | needs further information from LBL experiments. This is not an artifact of 
the expansion in Eq. fl5.27p . The reason is that and A^. T are small if \p s \ is small and 
this is true even if \q s \ is large, i.e. c s is large. This explains the limitations in our search 
for tests of the BBN constraint on c s derived in Section [Til]. On the other hand, in the 
present approximation, for the determination of the values of \p e \, \q e \ and \p s \ in principle 
the results of SBL experiments alone are sufficient. 

The free parameter \p e \ 2 is fixed by a measurement of the SBL v e oscillation 

amplitude. With the result of the LSND experiment and the negative result of all other SBL 
neutrino oscillation experiments |p e | 2 is approximately bounded by 0.5 x 10~ 3 < |j9 e | 2 < 10~ 2 . 
If the BBN constraint on c s is valid, then — > v T transitions are significantly suppressed in 
SBL neutrino oscillation experiments with an oscillation amplitude of order 10 -3 or smaller. 
Furthermore, Eq. (|5.30| ) shows that the oscillation amplitude A e . T is at most close to 0.1. 



In the discussion in Section [Til] we have seen that the BBN bound N v < 3.9 implies that 
c s <^ 10~ 2 . In this case we have 



6 Actually, this is always possible because of unitarity, but in general the relations are very 
complicated. 
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B ~ 



an d -<4/z;t — 7 ^./i;e ^e;r • 



(5.32) 



Inspecting Eq. (|5.27 ) we note that in scheme A the quantity 1— c M is now of the same order 
of magnitude as c e . As mentioned in the previous Section, this is particularly interesting for 
Amg BL < 0.3 eV 2 where the bound on B KfX from the disappearance experiments is not 
very stringent and, consequently, in this Am| BL region B^ should be much smaller than 
the present upper experimental upper bound. In addition, also 1 — c r is of the same order of 
magnitude, showing that at zeroth-order in the expansion over the small quantities c e and 
c s , the mixing matrix is given by 



/0 



u 



\ 



w 
w 








(5.33) 



Therefore, we have 



cos 9 z/3 + sin 9 



~ cos 7 V\ + sin 7 z/ 2 



— sin 9 V3 + cos 9 z/4 , 
t — sin 7 Vi + cos 7 v 2 



(5.34) 



Hence, the u e , v s - z/ 3 , z/ 4 and v^v T - v^-, v 2 sectors are decoupled and the oscillations of solar 
and atmospheric neutrinos are independent. Furthermore, since only the small mixing angle 
MSW solutions of the solar neutrino problem seems to be allowed by the data in the case 
of v e — > v s transitions (see Refs. [|3l],[^,[l(|]), at zeroth-order in the expansion over the small 
quantities c e , c s and 9 we have 



"3 



(5.35) 



Similar conclusions are valid in scheme B with c a replaced by 1 — c a and 1,2^3, 4. 

Let us now discuss LBL neutrino oscillations. Keeping only the dominant terms we get 



P 



(LBL) 



(-) C-) 



1 ~P, 



(LBL) 



(-) (-) 



sin 2 (2 7 ) sin 2 (0/2) 



(5.36) 



where ~ Am 2 1 L/2p in scheme A and Am 2 3 L/2p in scheme B. This shows that there must 
be significant LBL — > v T oscillations because of the sizeable disappearance known from 
the atmospheric neutrino anomaly, which implies that sin 2 (27) is large. For v e and 

v e — > v T transitions matter effects have to be taken into account. The upper bounds on 
these transition probabilities derived in Ref. [23 are of order 10 -2 . If c s is small, as follows 
from the BBN constraint, the "matter-stable" bounds on P, , , , and P^ LBL ] . derived in 



(-) (-) 



(-) (-) 



Ref. are improved for small values of Amg BL because of the inequalities 



C T > 1 - Ot 



leading to 



,(LBL) 

(-) (-) 



P^ L U < 2a°(l - a°) + c s (l - 2a°) 



(5.37) 



(5.38 
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Hence, the BBN constraint on c s imply that 



P (Tl) > P tTU <2ag(l-o2). (5.39) 

Obviously, the observation of a violation of these inequalities in LBL neutrino oscillation 
experiments would imply that the BBN constraint on c s is not valid. 



VI. CONCLUSIONS 

In this paper we have focused our discussion on the two four-neutrino schemes A and B 
(see Eq.( [1.5| )) which are compatible with the results of all neutrino oscillation experiments 
f22| , |23"|1 . As shown in Ref. ||23|| , the bound N v < 3.9 for the effective number N v of light 
neutrinos before the onset of BBN implies a stringent upper bound on c s = J2k=i2 \ U s k\ 2 
in scheme A and on 1 — c s in scheme B. Our limits are more stringent than those of Ref. 
23| because we used the complete collision rates presented in Ref. |]42|| . Thus, in scheme A 



we have obtained c s < 5 x 10 . We have demonstrated that in scheme B the limit derives 
from resonant transitions due to Am| BL , leading to an upper bound for 1 — c s that is more 
stringent than the one for c s in scheme A: 1 — c s < 5 x 10~ 5 . 

The validity of these limits from BBN depends to a certain extent on the results of 
the measurements of the primordial deuterium abundance which are controversial at the 



moment pyfl - We have therefore emphasized the importance of terrestrial experiments to 
get information on the value of c s in order to check the BBN constraints. In Section |TV] we 
have worked out lower bounds on c s in scheme A (1 — c s in scheme B) involving the following 
observables in SBL experiments: 

i. the v ^ — > v T oscillation amplitude A WT ; 

ii. the oscillation amplitude in disappearance experiments for Amg BL < 0.3 eV 2 ; 

iii. the oscillation amplitudes B e . e in reactor disappearance experiments and A^ e and A e . T 
in accelerator experiments. 

To proceed further, in Section [V] we have introduced a parameterization of the 4x4 



neutrino mixing matrix U (|5.7|) which is particularly suited for the schemes A and B. In this 
parameterization the quantity c s is given by c s = |p s | 2 + |g s | 2 in scheme A (1 — c s = \p s \ 2 + \<ls\ 2 
in scheme B), where p s and q s are in general complex parameters in U. The parameterization 



has lead us to the following observations: 

(a) The parameter \p s \ can be determined in SBL — > v e ,p,,r oscillation experiments. 

(b) For the determination of \q s \ in SBL oscillation experiments it is necessary to employ 
a u T neutrino beam and perform a disappearance experiment. 

(c) For N u < 3.9 the SBL oscillation amplitude A^. T is at most of order 10~ 3 whereas the 
LBL oscillation amplitude must be of order one in the same channel. 
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As already pointed out in Ref. ||23|| , the BBN constraint N v < 3.9 leads to a 4-neutrino 
mixing picture in which, apart from small corrections of the order of the small quantities c e 
and c s , the mixing matrix U is composed by the two decoupled sectors u e ,u s - z/ 3 ,z/ 4 and 
Up, v T - Ux, i>2 hi scheme A (y e , v s - u%, v 2 and v^, v T - u 3 , z/ 4 in scheme B). Consequently, the 
oscillations of solar and atmospheric neutrinos are independent. The solar neutrino problem 
is solved with u e — > v s oscillations governed by the mixing angle 9 and the atmospheric 
neutrino anomaly is explained with — > v T oscillations governed by the mixing angle 7 
(see Eq.( |5.7| ) for the definition of these angles). Concerning the solar neutrino deficit, it was 
shown that a vacuum oscillation solution involving sterile neutrinos is rather disfavored by 



the data (see Refs. [^T|,^]). Therefore, the specific prediction of a small mixing angle MSW 
solution of the solar neutrino problem with only v e — > v s transitions serves as an indirect 
check of the schemes A and B under consideration and of the BBN bound. This prediction 
implies a characteristic strong distortion of the 8 B v e spectrum and a day-night asymmetry 
in the neutral current event rate at SNO pH] . Experimentally, this prediction can be checked 



independently from solar model calculations (see Ref. |59| 



If indeed the BBN constraints on the 4-neutrino mixing matrix U are confirmed then the 
biggest gap in our knowledge of U, the sterile neutrino mixing, is considerably narrowed. In 
this paper we have shown that it is possible to test the ensuing mixing matrix to some extent 
in terrestrial experiments. The most striking feature in this context is the suppression of all 
SBL oscillation amplitudes. 
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FIGURES 
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FIG.0. Plot in the A 
the right of this curve (see Eq. ( |4.8| )) the lower bound ( |4.9| ) for c s applies, 
curve representing B® constitutes an upper bound for A^. T due to the conservation of 
probability. The values of a° and 5° have been obtained from the 90% CL exclusion 
curves of the CDHS |36| and CCFR |3^j experiments. The dash-dotted curve, dash-dot- 
dotted and dotted curves represent, respectively, the most recent exclusion curve presented 
by the CHORUS collaboration |52|] , the expected final sensitivity of the CHORUS |49[ and 
NOMAD [KO] experiments and the expected sensitivity of the COSMOS |ETJ experiment. 



FIG.0. Plot in the c s -Amg BL plane of the lower limit 2a° — 
the inequality ( [4.8|) . The values of a? e and 

curves of the Bugey |35|], CDHS pH] and CCFR p7[ experiments. 
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FIG.[3[ Plot in the A^. T -/S.m\ Bh plane of the upper bound (a°) 2 (solid curve) for A 
Eq.( 4.11|) ) in the case of a very small c s . The values of a® have been obtained from the 90% 
CL exclusion curve of the Bugey experiment. The dash-dotted curve, dash-dot-dotted 
and dotted curves are the same as in Fig.[TJ 
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